mode control of quantum systems 



Daoyi Dong 

School of Information Technology and Electrical Engineering, University of New 
South Wales at the Australian Defence Force Academy, Canberra, ACT 2600, 
Australia. 

E-mail: daoyidong@gmail.com 

Ian R. Petersen 

School of Information Technology and Electrical Engineering, University of New 
South Wales at the Australian Defence Force Academy, Canberra, ACT 2600, 
Australia. 

E-mail: i.r.petersen@gmail.com 

Abstract. This paper proposes a new robust control method for quantum systems 
with uncertainties involving sliding mode control (SMC). Sliding mode control is a 
widely used approach in classical control theory and industrial applications. We show 
that SMC is also a useful method for robust control of quantum systems. In this paper, 
we define two specific classes of sliding modes (i.e., eigenstates and state subspaces) 
and propose two novel methods combining unitary control and periodic projective 
measurements for the design of quantum sliding mode control systems. Two examples 
including a two-level system and a three-level system are presented to demonstrate 
the proposed SMC method. One of main features of the proposed method is that 
the designed control laws can guarantee desired control performance in the presence 
of uncertainties in the system Hamiltonian. This sliding mode control approach 
provides a useful control theoretic tool for robust quantum information processing 
with uncertainties. 
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1. Introduction 

Controlling quantum phenomena has been an implicit goal even since the establishment 
of quantum mechanics |l]-[3]. Many practical tasks arising from atomic physics [3], 
molecular chemistry [2], [l]-[7] and quantum optics [8] can be formulated as quantum 
control problems. It is desired to develop quantum control theory to establish a firm 
theoretic footing for the active control of quantum systems. Quantum control theory 
has played an important role in guiding the control of chemical reactions p], [7]. 
Recently, the development of quantum control theory has been recognized as a key 
task required to establish practical quantum information technology [8], [9], [10], [TT] . 
Some useful tools from classical control theory, such as optimal control theory [12] , [13] 
and feedback control approaches [Il]-[22] have been applied to problems of population 
transfer, quantum state preparation, quantum error correction and decoherence control 

Although some useful results have already been obtained, research on quantum 
control is still in its infancy. From the perspective of practical applications, it is 
inevitable that there exist uncertainties and disturbances in external fields, system 
Hamiltonians or initial states [2l|. Many cases of unknown information and errors can 
also be treated as uncertainties. Hence, the requirement of a certain degree of robustness 
in the presence of uncertainties and noises has been recognized as one of the key aspects 
for developing practical quantum technologies [25]- [35]. Several approaches have been 
introduced to enhance the robustness of quantum systems. For example, Zhang and 
Rabitz [26] used a minmax approach to analyze the robustness of molecular systems. 
In [36], James and coworkers have formulated and solved a quantum robust control 
problem using a method for linear quantum stochastic systems. In this paper, we 
develop a new sliding mode control approach to enhance the robustness of quantum 
systems. Sliding mode control (SMC) is a useful robust control strategy in classical 
control theory and industrial applications [37|-[39|. References [IQ] and [H] have briefly 
discussed the possible application of sliding mode control to quantum systems. This 
paper will formally present sliding mode control methods for quantum systems to deal 
with uncertainties. 

The sliding mode control approach generally includes two main steps: selecting a 
sliding surface (sliding mode) and controlling the system state to and maintaining it in 
this sliding surface. The sliding surface guarantees that the quantum system has desired 
dynamics in this surface. We will select an eigenstate of the free Hamiltonian or a state 
subspace of the controlled quantum system as a sliding mode. To control the system 
state to and then maintain it in this sliding surface, most feedback control methods 
in classical control theory are not directly applicable since we generally cannot acquire 
measurement feedback information without destroying the quantum system state. It is 
necessary to develop new approaches to accomplish this task. Hence, we propose two 
new unique methods for this task: one method involves combining time-optimal control 
design and periodic projective measurements for the case where an eigenstate is taken 
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as the sliding mode; the other method is implemented through quantum amplitude 
amplification and periodic projective measurements for the case where a state subspace 
is identified as a sliding mode. It is worth noting that an important assumption required 
in the proposed methods is that accurate projective measurements (using the eigenstates 
of the free Hamiltonian as the measurement basis) are possible. A connection to the 
quantum Zeno effect is briefly discussed and two specific examples are presented to 
demonstrate the proposed methods. The main feature of our methods is their robustness 
to uncertainties in system Hamiltonians. The sliding mode control method has potential 
applications to preparation of quantum states and quantum error correction. 

This paper is organized as follows. Section 2 defines two classes of candidate sliding 
modes and formulates the two quantum control problems considered in this paper. In 
Section 3 we present a sliding mode control method combining time-optimal control 
design and periodic projective measurements for quantum systems where an eigenstate 
is identified as a sliding mode. An illustrative example of a two-level quantum system 
is analyzed in detail. Section 4 proposes a control design method using quantum 
amplitude amplification and periodic projective measurements for controlling quantum 
systems when the sliding mode is a state subspace. An example of a three-level system 
is presented to demonstrate the proposed method. Concluding remarks are given in 
Section 5. 

2. Sliding modes and problem formulation 

In this paper, we consider finite dimensional (A^-level) quantum systems whose dynamic 
evolution is described as follows (setting h = I'm. this paper) 

i^^\m) = Ho\m) + Y.ukmk\m); m = o)) = i^o), (i) 

k 

where Hq is the free Hamiltonian, Uk{t) G R is the external control, and {Hk} is a 
set of time-independent Hamiltonian operators. If the eigenvalues and corresponding 
eigenstates of Hq are denoted as \j and \4>j){j = 1,2, . . . , N) (i.e., -f^ol^j) = 
respectively, \ipit)) can be expanded as: 

N 

\m) = J2^jm,). (2) 

i=i 

In ([1]), H = Hq + J2k'^k{t)Hk is Hermitian which ensures that the transition 
from a pure state |'?/'(0)) to another I'lpit)) can be accomplished through a unitary 
transformation U{t)] i.e., \ip{t)) = U{t)\ip{0)). The control problem is converted into 
the problem that given an initial state and a target state, find a set of controls {uk{t)} 
in ([T]) or a unitary transformation U{t) to drive the controlled system from the initial 
state into the target state. 
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2.1. Sliding modes 

Sliding mode control (SMC) is a useful approach to robust controller design for 
electromechanical systems [39]. In SMC, a sliding surface (sliding mode) and a switching 
function are determined [38]. A sliding mode is generally defined as a specified state 
region where the system has desired dynamic behavior. A switching function is designed 
to control the system state to the sliding mode. SMC has two main advantages |38j : 
(i) the dynamic behavior of the system may be determined by the particular choice 
of switching function; (ii) the closed-loop response becomes totally insensitive to a 
particular class of uncertainties. It is the order reduction property and low sensitivity to 
uncertainty that makes SMC an efficient tool for controlling complex high-order systems 
subject to uncertainties [39] . 

To apply the idea of SMC to quantum systems, we first need to define a sliding mode 
where the quantum system has desired dynamics. A sliding mode can be represented 
as a functional of the state and Hamiltonian H\ i.e., S{\il)),H) = 0. For example, 
an eigenstate {(pj) of Hq can be selected as a sliding surface. In this case, we can 
define S{\^p),H) = 1 - If the initial state |^o) 

is in the sliding mode; i.e., 

S{\4'o),H) = 1 — \{i'o\<Pj)\'^ = 0, we can easily prove that the quantum system will 
maintain its state in this surface under the action of the free Hamiltonian Hq only. In 
fact, \^p{t)) = e-^-^«*|^o) and 

s{\m),H) = 1 - mm)\' = i - me'^'^'m' = i - k^oI^.^^^^t 

= l-|(^o|0,)ne^'^f = 0. 

That is, an eigenstate of Hq can be identified as a sliding mode. 

More generally, an invariant state subspace of a quantum system satisfying 
S{\tp),H) = can be defined as a sliding mode. For example, the wavefunction 
controllable subspace considered in [ID] can be identified as a sliding mode. Considering 
a simple control u{t) (i.e.. A; = 1 in ([1])) and substituting ([2]) into ([Tj), we denote 
C{t) = {ci(t), C2(t), . . . , C7v(t)} and can obtain [12], [13]: 

iC{t) = AC{t) + u{t)BC{t); C{t = 0) = Co; (3) 

N 

Co = (coj)jli; CQj = (0j#o), I^ojI^ = ^' 

where A and B correspond to the operators Hq and Hi (i.e., /c = 1 in ([T])), respectively. 
Consider the Model I presented in [ID], 

iC{t) = AC{t) + u{t)BC{t) = {A + u{t)B)C{ty, C{t = 0) = Co, (5) 

where the free Hamiltonian of the five-level system is A = diagjl.0, 1.2, 1.3, 2.0, 2.15} 
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and the control Hamiltonian H,, 
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u(t)B is as follows: 

uit)\ 



uit) 

J 



(6) 



In [lO], we have proven that the subspace spanned by {\4>i), |04), l^s)} is a 
wavefunction controllable subspace [15]. We may select Q as a. sliding surface. We 
can easily prove that if the initial state of this system is in this sliding surface, its state 
will be maintained in this surface under the action of Hamiltonian H = A + u{t)B. In 
fact, we can express the sliding mode as follows: 
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After straightforward calculations we obtain the following relationship : 
d 



dt 



(|ci(t)|^+|c4(t)|^ + |c5(t)h 



0. 



it is clear that 
ation from ([5]): 

(7) 



(8) 



From this, it is clear that |ci(i(:)p + |c4(t)p + |c5(t)p = |coip + |co4p + |co5p = 1. That is, 
S{\il!{t)), H) = 0; i.e., the state of this system will be maintained in this surface under 
the action of Hamiltonian H = A + u{t)B. 

Some other state subspaces such as decoherence-free subspace jlB], [IZ] can also be 
defined as sliding modes. However, this paper will focus on the two particular classes 
of sliding modes. The first class corresponds to a sliding mode which is an eigenstate. 
The second class corresponds to a sliding mode which is a wavefunction controllable 
subspace. 



2.2. Problem formulation 

In the above subsection, we have presented two candidate sliding modes. If a quantum 
system state is driven into a sliding mode, the state will be maintained in the sliding 
surface under the action of some class of Hamiltonians determined by the sliding 
mode. However, in practical applications, it is inevitable that there exist noises and 
uncertainties in system Hamiltonians, initial states or control fields. An important 
advantage of sliding mode control is its robustness against uncertainties. Our main 
motivation of introducing sliding mode control to quantum systems is to deal with 
these uncertainties. In this paper, we focus on the class of uncertainties which can be 
approximated as a perturbation in the system Hamiltonian. For example, the unitary 
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error in and the operational error in a quantum logic gate can be classified into this 
class of uncertainties. We also suppose that the uncertainties are bounded. We represent 
the uncertainties as i^A = Yl:i^i{^)Hi-, where ei{t) G R, a/^T^TW < ^ (e G R"^) and 
{Hi] is a set of time-independent Hamiltonian operators. We further suppose that the 
system is completely controllable [IS]. 

The control problem under consideration is stated as follows: for a given initial 
state, design a control law to steer the quantum system state into and then maintain 
the state in a sliding mode domain (a neighborhood containing the sliding mode) in the 
presence of bounded uncertainties in the system Hamiltonian. This quantum sliding 
mode control problem is greatly different from the traditional sliding mode control. 
Once the uncertainties take the state slightly away from the sliding mode, there is 
always a finite probability (we call it the probability of failure) that the system state 
will collapse out of the sliding mode domain when one makes a measurement on this 
system. Hence, if the allowed probability of failure is we may define the sliding mode 
domain V = : KV^|$)P > 1 -po, l-^*) e {|$) : S{\^),H) = 0}}. We expect that 
the control law can ensure that the system state remains in the sliding mode domain P 
except that a measurement operation may take it away from V with a small probability 
(not greater than po). The quantum sliding mode control problem includes three main 
subtasks: (I) for any initial state (assumed to be known), design a control law to drive 
the system state into a defined sliding mode domain V; (II) design a control law to 
maintain the system state in T>] (HI) design a control law to drive the system state 
back to P if a measurement takes it away from V. For convenience, we suppose that 
there exist no uncertainties during the control processes (I) and (HI). In particular, we 
consider the following two quantum control problems. 

(1) Quantum Control Problem 1 (QCPl): For an uncertain quantum system in 
which an eigenstate \(f)j) of Hq defines a sliding mode (i.e., S{\ip), H) = 1 — \ {ip{t)\(f)j)\'^ = 
0), (I) drive the controlled quantum system state into the sliding mode domain 
^ = ili^) '■ — 1 ~ Po}j (II) maintain the system state in V except that 
measurements may take it away from V with at most probability po] (III) if the system 
state is taken away from V, design a control law to drive it back to V. 

(2) Quantum Control Problem 2 (QCP2): For an uncertain quantum system in 
which a wavefunction controllable subspace Q defines a sliding mode (i.e., S{\ip),H) = 
1 — |(-?/^(t)|$)p = where 1$) G Q), (I) drive the controlled quantum system state into 
the sliding mode domain V = {lip) : K'?/'!^)^ > 1 —po, |$) G Q}; (II) maintain the state 
in V except that measurements may take it away from T> with at most probability po] 
(III) if the system state is taken away from V, design a control law to drive it back to 
V. 
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3. SMC based on time-optimal design and periodic measurements 

3.1. The general method 

In this section, we consider QCPl, in which an eigenstate is identified as a sliding mode. 
Our first task is to design a control law to drive the system state to this chosen sliding 
surface. Since we ignore the effect of uncertainties during the control processes (1) and 
(111), we wish to accomplish this task as quickly as possible. Hence, we will use a time- 
optimal control approach for this task. For the subtask (II), we use periodic projective 
measurements to achieve our goal. In coherent control, measurement is usually regarded 
as having deleterious effects. Recent results have shown that quantum measurements 
can be combined with unitary transformations to complete some quantum manipulation 
tasks and enhance the capability of quantum control [H], [19], [50]- [55]. For example, 
Vilela Mendes and Man'ko [H] showed that nonunitarily controllable systems might 
become controllable by using "measurement plus evolution" . Roa et al. [53] have used 
sequential measurements to control quantum systems. Rabitz and coworkers [50]-[52] 
have demonstrated that projective measurements can serve as a control tool. In this 
section, we will combine time-optimal control and periodic projective measurements to 
accomplish sliding mode control of quantum systems. 

The steps of the control algorithm for QCPl follows: 

(i) Select an eigenstate of Hq as a shding mode S{\ip),H) = 0; 

(ii) For a known initial state \ipo), design a time-optimal control law that can drive \ipo) 
to the sliding mode S; 

(iii) For eigenstates \(f)k) {k 7^ j), design corresponding time-optimal control laws that 
can drive to S using a similar method to that in (ii); 

(iv) For given po and e, design the period T for the projective measurements; 

(v) Use the designed control law to drive the system state to 5*, then implement 
periodic projective measurements with the period T to maintain the system state in 
V = {\tp) : KV'l^i)!^ > 1 ~Po}- If the state collapses to 10^) due to a measurement, 
we use the corresponding control law to drive it to S and then continue to make 
periodic projective measurements. 

From the above, we can see that the design of a time-optimal control law and the 
selection of the period T for projective measurements are the two most important tasks 
in this control algorithm. The time-optimal control problem has been an interesting 
topic in quantum control, in which it is required to design a control law to achieve a 
desired state transfer in a minimum time in order to minimize the effects of relaxation 
and decoherence [13], [56]. Khaneja and coworkers [13] have studied time-optimal control 
of spin systems under the assumption of unbounded controls. Boscain and Mason [57] 
have investigated the time-optimal control problem for a spin 1/2 system with bounded 
controls. For several simple quantum systems (e.g., two-level systems), it is possible to 
obtain analytical results. However, it is generally difficult to find a complete solution 
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for high-dimensional quantum systems. In these cases, it may be useful to develop a 
numerical simulation method to find an approximate solution. 

Another important task is to design the measurement period T so that the control 
law can guarantee control performance. An extreme case occurs when T — 0. 
That is, after the quantum system state is driven into the sliding mode, we make 
frequent measurements. This corresponds to the quantum Zeno effect [58] which can 
guarantee that the state is maintained in the sliding mode in spite of the existence 
of uncertainties. However, it is a difficult task to make such frequent measurements 
in practical applications. We may think that the smaller T is, the bigger the cost 
of measurements becomes. Hence, we wish to design a period T which is as large as 
possible, when we have a bound e on the uncertainties and require a probability of 
failure po- In the following subsection, we will present a specific example of a two- level 
quantum system to demonstrate how to design the period T. 

Remark 1 It is clear that the sliding mode control approach can he used in the 
preparation and protection of quantum states under uncertainty conditions. Moreover, 
if the initial state is unknown, we can first make a projective measurement. If the result 
is the eigenstate \(t>j), we would continue to implement periodic projective measurements 
to maintain the system state in V. Otherwise, we would use a corresponding control 
law to drive the system state to S. This slight amendment to our approach enables our 
method to achieve robustness against variations in the initial state as well as robustness 
to uncertainties in the system Hamiltonian. 



3.2. An illustrative example: two-level system 

To demonstrate the proposed method, here we consider a two-level quantum system, 
which can be used as a quantum bit (qubit) and has important potential applications in 
quantum information. In practical applications, we often use the density operator p to 
describe the state of a quantum system. For a pure state the corresponding density 
operator is p = \ip){ip\. For a two-level quantum system, the state p can be represented 
in terms of the Bloch vector r = {x,y,z) = (trlpcr^.}, tr{po"j^}, trjpcr^}): 

p = l(/ + r-a). (9) 

where a = {ax, o-y,az) are the Pauli matrices described as follows: 

The evolution of p can be described by the equation 

p = -z[H,p], (11) 

where H = Hq + Y.k=x,y,z^k{t)Ik, Hq = 1^ = ^a^, 4 = and ly = \(yy. 

Without loss of generality, we select the sliding mode as S'(| ■?/'), -f^) = 1 — I ("^10)1^ = 0. 
This means that we select the eigenstate |0) of Hq as the sliding mode. If we have 
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driven the system state to the shding mode at time to 5 it will be maintained in 
this sliding mode using only the free Hamiltonian Hq, i.e., S{\ip(t>t(,)) ^ Hq) = 0. We 
assume that the possible uncertainties in the system Hamiltonian are represented by 
Ha = e^{t)h + e,{t)h, where el{t) + el{t) < e. Now we use the control algorithm in 
Section 3.1 to accomplish the robust control design. 

For simplicity, we assume \ipo) = |1) (this makes the subtasks (I) and (HI) in QCPl 
become the same) and consider = u{t)Ir^. If u{t) is not bounded, it is convenient to 
design a time-optimal control law to drive |1) to the sliding mode |0) using the results 
in [13], [59]. From Theorem 1 in [59J, we learn that the minimum time required to 
accomplish this task is 0. If u{t) is bounded (i.e., \u{t)\ <V,V& R-'*'), we can use the 
method in [57] to design a time-optimal control law to drive |1) to |0). 

Now we consider the design of the measurement period T. First consider an 
uncertainty represented by = ez(t)Iz (where \ez(t)\ < e). If S{\ipo),H) = 0, for 
Ha = ^zh we have 

S{\m,H) = 1 - Hmm' = l - |(^o|e^(^«+^^'^)*|0)|^ = l - |(^o|0)r|e^(5+f-)f = o. 

This uncertainty does not take the system state away from the sliding mode. Hence, we 
may ignore this uncertainty in the following analysis. 

Now consider a bit-flip type uncertainty of the form e(t)Ix ( |e(t)| < e). We have 
the following theorem. 

Theorem 1 For a two-level quantum system with the initial state (xq, 1/0,-^0) = 
(0,0,1) (i.e., \0)), the system evolves to (x", j/f, 2;") and {x\,y\,z\) under the action of 
the Hamiltonians if" = Iz + e(t)Ix (where \e{t)\ < e) and = Iz + llx (where e = +e 
or —e), respectively. Then for arbitrary t G [0, ^^^-^ ] ; > z\. 



The proof of this theorem is presented in Appendix. Theorem 1 shows that z\ 
can be taken as an estimate of the bound on z^. From the proof, it follows that when 

—2 — — — . 

t = ^q^ y the probability of failure is p' = If Po < p', using fl51|) in Appendix we 
can choose T according to the following relationship: 



l-ZT 1- cos VT+i^T 

P' = ^ = 2 • (^2) 

Hence, we may choose the measurement period T as follows: 

T = ) 7^ arccos[l - if Po < ^] ^^3^ 



otherwise. 



In quantum computation, an important result is the fact that arbitrarily accurate 
quantum computation is possible provided that the error per operation is below a 
threshold value [60]. In this application the uncertainties under consideration may 
come from quantum gate errors. If we define the gate fidelity as follows [9], [30] : 

r{u,,u) = mmmuium, 
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a straightforward calculation shows that the gate fidelity is not less than -^7^=7 under 
our control strategy. If we define the quantum gate error as = 1 — J-'{Uo,U), the 
gate error is not greater than 1 — -^tj^- Hence, the proposed method can be used in 
the design of robust quantum gates. 

Now we consider a specific case. Suppose that the control is bounded \u{t)\ < 10, 
the bound of uncertainties is e = 0.1, and the allowed probability of failure is po = 1.00%. 
According to [57], the time-optimal control from |1) to |0) is bang-bang control and the 
number of switchings required is 1. A straightforward calculation using the method 
in [57] leads us to the conclusion that we use u(t) = —10 in t G [0,0.1573] and use 

—2 

u{t) = 10 in t G (0.1573,0.3146]. Since po > y+F' using ( fT3l) we can obtain the 
measurement period T = 3.1260. The corresponding maximum gate error is not greater 
than Ge = l- ^ = 0.50%. 

4. SMC based on amplitude amplification and periodic measurements 

In this section, we identify a state subspace as a sliding mode and employ the quantum 
amplitude amplification method to design the control laws. We will first introduce 
the quantum amplitude amplification method, then present the control algorithm, and 
finally give an illustrative example. 



4.I. Amplitude amplification 

The quantum amplitude amplification method is a powerful approach used in many 
quantum algorithms [61]- [M]. The central task in quantum amplitude amplification is 
to find a suitable operator Q whose repeated action on the initial state can increase 
the probability of chosen eigenstates. If we denote X = {|0), . . . , . . . , — 1)} as a 
set of orthonormal basis in the A^-dimensional complex Hilbert space Ti, a pure state 
lip) of an A^-level quantum system can be represented as \ip) = J2x=o^x\x), where 
'^x=o l*^^'!^ = 1. A Boolean function x : X ^ {0, 1} defines two orthogonal subspaces of 
H: the "good" subspace and the "bad" subspace. The good subspace is spanned by the 
set of basis states \x) G X satisfying = 1 and the bad subspace is its orthogonal 
complement in H. We may decompose \ip) as lip) = \ipg) + \ipb), where \ipg) = Pg\ip) 
denotes the projection of \ip) onto the good subspace with the corresponding projector 
Pg, and lipb) = {I ~ Pg)\i^) denotes the projection of onto the bad subspace (here 
/ is the identity matrix). It is clear that the occurrence probabilities of a "good" state 
\x) [x(yX) = 1] and a "bad" state \x) [x(x) = 0] upon measuring \ip) are g = {'ipglipg) and 
b = {ipbl'ipb) = ^ — 9, respectively. 

Let lip) = U\0). Given two angles < (pi, < tt, quantum amplitude amplification 
can be realized by the following operator [61] 

Q = Q(W,x,<^i,<^2) = -UV^'U-'V^\ (14) 
The operators Vq^ and V^^ conditionally change the phases of state |0) and the good 
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states, respectively |61j, and they can be expressed as [52] : 

Po^^=/-(l-e^^^)|0)(0|, (15) 
=/- \x){x\. (16) 

X(x)=l 

The action of Q can be described by the following relationship |52j : 
Q|^)= [(i_e*c,i)(i_^_^^g.^.)_g.^.]|^^^ 

+ [^(1 -e*^i)(e*^' - 1) -e^^'ll^b)- 
Thus, we can amplify (or shrink) the amplitude of \iljg) (or \ipb)) by a suitable selection 
of the parameters (^2 in Q. 



4-.2. The control algorithm 

The main steps in sliding mode control based on amplitude amplification and periodic 
projective measurements for QCP2 are as follows: 

(i) Select a state subspace f2 as a sliding mode S{\ip), H) = 0; 

(ii) For a known initial state \ipo)y identify as a "good" subspace and construct an 
amplitude amplification operator Qd^/'o), S) to amplify the probability of projecting 
iV'o) into fl; 

(iii) Using the probability po and CI{\iIjo),S), determine a number Lq of iterations 
required to guarantee that the control law drives the system state into the sliding 
mode domain V = : KV'I^)^ > I - Po, \^) e 

(iv) For other eigenstates \(f)k) not in the "good" subspace, first apply a unitary 
transformation Uk on \(f)k) to obtain a superposition state {ipk)- Then construct a 
corresponding amplitude amplification operator Q{\ipk), S) and choose the required 
number Lk of iterations using a similar method as used in (ii) and (iii); 

(v) Using pq and e, choose the period T for the periodic projective measurements; 

(vi) Use the designed control law to drive the system state into V, then implement 
periodic projective measurements with the period T to maintain the system state 
in V. If a measurement result is \(j)k) which is not in Q, use the corresponding 
control law in (iv) to drive the system state into T>. 

Remark 2 Here we use amplitude amplification as an important part of our control 
algorithm for the sliding mode control of quantum systems. This is similar to the idea 
in 15^ . However, 152^ does not consider the issue of robustness. Here, our goal is 
to develop a new method to deal with uncertainties in the system Hamiltonian. From 
the above control algorithm, it is clear that we can design different controllers offline. 
Hence, this is a convenient approach to be applied to different control tasks. 

Remark 3 In the above control algorithm, the construction of amplitude amplification 
operator Q is dependent on the initial state; i.e., the initial state must be known. When 
the initial state is an eigenstate \4>k) , we. usually need apply a unitary transformation 
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Uk to drive to a superposition state before constructing Q. Here Uk may be a small 
perturbation or an easily-realized unitary transformation. If the initial state is unknown 
or is a mixed state, we can construct a Kraus-map to control the initial state to a specified 
pure state (for details, see, e.g., l65^ ) and then use the proposed control algorithm to 
accomplish the control task. 

Remark 4 In quantum information, quantum error correction is an important problem. 
Two important paradigms for quantum error correction have been proposed: active error 
correction and passive error avoidance. The main idea for active error correction is 
encoding quantum information using redundant qubits, detecting possible errors and 
then correcting the error JEM/, W^ , In the paradigm of error avoidance, one may 
encode quantum information in a decoherence-free (noiseless) subspace J2^ , fj^ , W^ - 
Our control method presented above provides a possible unified framework for the two 
paradigms of quantum error correction. We may identify a decoherence-free subspace 
as a sliding mode. In the sliding mode, the system is robust against some specified 
uncertainties (errors). However, some other classes of errors may break the specific 
symmetry in the system- environment interaction and take the system state away from 
the decoherence-free subspace. In that case, we can design a control law to drive the 
system state back to the sliding mode. A schematic demonstration of the relationships 
between sliding mode control (SMC) and quantum error correction (QEC) is shown in 
Figure [I]. From this discussion it can be seen that SMC has the potential to provide a 
new method for quantum error correction. The detailed development of this idea is an 
area for future research. 



4.3. An illustrative example: three-level system 



We now consider a three-level quantum system. In ([3]), let u{t) = 1 and A 
diag{ — 1, 0, 1}. The matrix B and the uncertainty matrix follows: 

/ 1 \ / e(t) \ 



B 



1 





















18) 



where |e(t)| < e. This leads to the following equation of the form ([3]) 

/ ci(t) \ /-II <t) \ I ci{t) \ 
i C2{t) =10 C2{t) . (19) 

V csit) J \ e{t) 1 J \ c,{t) J 

The three eigenstates {|0i), \(p2), l^s)} are denoted as {|1), |2), |3)}, respectively. 
Suppose that the initial state is |1); i.e., {ci(0), 02(0), 03(0)} = {1,0,0}. The three- 
level system corresponding to (|T8l) can be represented schematically as shown in Figure 
[21 The form of the B matrix implies that there is a direct coupling between states 
|1) and 1 2). The form of the uncertainty matrix implies that a disturbance 
coupling e(t) exists between |1) and |3). It is easy to check using the method in 



Dong and Petersen, SMC for Focus on Quantum Control 



13 



r" 



SMC 



33- CL 




I 

Uncertain^ ^C - f — !> Error 




QEC 



CE < ; 



Figure 1. A representation of the relationships between sliding mode control (SMC) 
in QCP2 and a unified framework of quantum error correction (QEC). The sliding 
mode (SM Hsm) corresponds to the decoherence-free subspace (DPS); the complement 
space Hgj^i corresponds to the effect of the environment in QEC; and the uncertainties 
in SMC correspond to possible errors in QEC. The subtask (I) in SMC corresponds 
to the initialization of QEC; the subtask (II) in SMC corresponds to the protection of 
DFS in QEC; and the control law (CL) in the subtask (III) of SMC corresponds to the 
process of correcting errors (CE) in QEC. 




Figure 2. A three-level system. The subspace spanned by {|1),|2)} corresponds to 
the sliding mode (SM). The uncertainty e{t) exists between |1) and |3). 



Section 2.1 that the subspace spanned by {|1), |2)} can be used as a shding mode; i.e., 
S{\%Ij),H) = 1 — (K'i/'ll)!^ + |('?/'|2)p) = 0. Hence, our control problem may correspond 
to a practical problem in quantum information. For example, the two-level system with 
levels |1) and |2) can be used a qubit, and the uncertainty matrix H/\ may represent a 
possible leakage into states outside the qubit subspace [68] . 

We now consider the design of period T for the periodic projective measurements. 
Let Cj{t) = Xj(t) +iyj{t){j = 1,2,3), where Xj(t),yj{t) e R. Also, let the vector r]{t) 
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denote ■r]{t) = {xi{t),yi{t),X2{t),y2{t),x^{t),y?Xt)Y ■ From we obtain 
m = F{t{t))^{t) 



(20) 



where 



mt)) 



\ 






-1 





1 





1 





-1 





-e(t) 





1 










-1 

















e{t) 











-e{t) 











-1 







1 





(21) 



and 7^(0) = (a;i(0), yi(0), X2(0), 1/2(0), ^3(0), 2/3(0))^ = (1, 0, 0, 0, 0, 0)'^. 

We now introduce the co-state vector \{t) = (Ai(t), X2(t), Xsit), X^it), X^it), Xe(t))'^, 
and obtain the corresponding Hamiltonian function as follows: 

m{rj{t), e{t), A(t), t) = X'^{t)F{e{tMt) = e{t)M{t) + N{t), (22) 

where M{t) = Xr,{t)yi{t) - Xeit)xiit) - X2it)xsit) + Ai(t)y3(t) and N{t) = A2(t)xi(t) - 
Xi{t)xi{t) - Xi{t)yi{t) + X3{t)yi{t) - X2{t)x2{t) + Xi{t)y2{t) - Xeit)x3{t) + XrMvsit). 
According to Pontryagin's minimum principle [69], a necessary condition for e{t) to 
minimize the functional J(e) = xl{tf) + y^itf) is 

H(77(t),e(t),A(t),t) <e(77(t),e(t),A(t),t). (23) 

Hence, 

l(t) = -esgnM(t). (24) 

That is, the optimal control is a bang-bang control strategy; i.e., e(t) = e = +e or — e. 
Without loss of generality, now we let e{t) = e and focus on 

77(t) = F(e)r/(t), (25) 

where r/(0) = (1,0,0,0,0,0)^. 

Consider the optimal control with a fixed final time tf and a free final state 
Vitf) = {x,{tf),y,{tf),X2{tf),y2{tf),x,{tf),y,{tf)f. Let J{t) = xl{t) + yl{t). 
According to Pontryagin's minimum principle, A(t/) = ^J(t)\t=tf From this, we 
obtain (Ai(ty),A2(t/),A3(t/),A4(t/),A5(t/),A6(t/)) = (0,0,0,0,2a;3(t/),22/3(t^)). Now 
we consider another necessary condition X{t) = 
following relationships: 

X{t) = F{e)X{t). (26) 

By simulation of fl25|) and fl26|) . for different values of e, we find the sign of M{t) does 
not change over the first interval for which J{t) is monotonically increasing with respect 
to t. Hence we can estimate the required period T by considering the first interval on 
which J(t) is monotonic. For example, if we consider e = 0.1, we obtain the simulation 
results shown in Figure [3l From this, we obtain tf = 1.1160, X3{tf) = —0.0055, 
y^(tf) = -0.0705 and J{tf) = 0.0050. Then we can check that M(t) indeed does 



or] 
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t 



Figure 3. The curve of J(t) = x\{t) + y1{t) corresponding to ([25|) . 



0.06 




t 



Figure 4. The curve of M{t) corresponding to the first monotonic interval of J{t) in 
Figure O 



not change sign on t G [0,t/] as shown in Figure HI Now we can design the period T 
using the following relationship: 

^ ^ I ti (where ti G [0,t/) and J(ti) =po) iipo<J{tf)\ 
It/ otherwise. 

For example, if pq = 0.50%, we have po < J(tf). From Figure [3] and (p71) . we obtain 
T = 1.1160. 

Now we design a control law for the subtask (III) using the quantum amplitude 
amplification method introduced in Section 4.1. For simplicity, we use a simple choice 
of Q where ipi = n, (p2 = t'' and the subspace spanned by {|1), |2)} corresponds 
to the "good" subspace. When the measurement result is \(j)k) = |3), we first 
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apply a unitary transformation (a perturbation) W3 on |3) to obtain j-^s) = U^\'i) = 
0.0600|1) + 0.0800|2) + 0.9950|3). Now, using a similar method to the example in |52], 
we obtain the state |^^) = 0.5986|1) + 0.7981 12) + 0.0682|3) after 7 iterations of Q. Now 
we make a measurement on l^/^g), and the probability of failure is j9q = 0.47% (< 0.50%). 

5. Concluding remarks 

This paper focuses on a robust control problem for quantum systems with uncertainties. 
The main contributions are as follows: (1) We propose a new robust control framework 
based on sliding mode control to deal with uncertainties in quantum control. Sliding 
mode control is a powerful methodology in classical control theory and has many 
industrial apphcations. Hence, the proposed method has the potential to open a new 
avenue for robust control design of quantum systems, which will provide a new useful 
tool for quantum information processing with uncertainties. (2) This paper presents two 
specific control algorithms for two classes of quantum control problems. In particular, 
several approaches such as time-optimal control, projective measurements and quantum 
amplitude amplification are combined to design quantum control algorithms. Two 
examples are also analyzed in detail using the proposed algorithms. 

It should be pointed out that the results in this paper are only a first step in 
developing systematic sliding mode control approaches for quantum systems. Much 
further work is required and several open problems are listed as follows: (1) The 
measurement period T designed in this paper is relatively conservative and it may be 
possible to obtain a larger value of T when po is relatively large. (2) The connection 
between the proposed methods and quantum error correction and quantum state 
preparation has been discussed. However, the specific applications still require to be 
developed in detail. Moreover, we use two model systems as the examples to illustrate 
the proposed methods. A future task is to connect such model systems to specific 
quantum systems (e.g., quantum optical systems and spin systems) and to consider the 
experimental implementation of the proposed methods for such real quantum systems. 
(3) We only consider a simple amplitude amplification operator Q in our example, and 
the questions of how to determine an optimal Q and how to realize Q using specific 
control fields are still open problems. (4) We have ignored the uncertainties during the 
control processes (I) and (HI), which require a relatively short time to accomplish the 
two control objectives. An improvement worth exploring is to consider the effect of 
uncertainties in all three subtasks. (5) We have constrained the type of uncertainties to 
be dealt with, so future research could be directed forwards extending the sliding mode 
control approaches to deal with other types of uncertainties for quantum systems. 
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Appendix: The proof of Theorem 1 

For H = Iz + e(t)Ix, using p = —i[H, p] and Qj, we have 
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(28) 



We now consider e{t) as a control input and select the performance measure as 

J{e) = zt,. (29) 

Also, we introduce the co-state vector \{t) = {Xi{t) , X2it) , X^it))'^ and obtain the 
corresponding Hamiltonian function as follows [69j: 

e(r(t),e(t),A(t),t) = -Xi{t)yt + X2{t)xt + e{t){X3{t)yt - X2{t)zt). (30) 

According to Pontryagin's minimum principle [69], a necessary condition for i{t) to 
minimize J(e) is 

H(f (t), e(t), X{t),t) < e(f (t), e{t), X{t),t). (31) 

Hence the optimal control i{t) should be chosen as follows: 

e(t) = -esgn(A3(t)yt - X2{t)zt). (32) 

That is, the optimal control strategy for e{t) is bang-bang control; i.e., e{t) = e = 
+e or — e. Now we consider the system 

-1 







yt 
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(33) 



i+e 



cos tut + 
sin cut 



-2 ^ 1 

" ■ COS out + 



(34) 



Now consider the optimal control problem with a fixed final time tf and a free 
final state r^^ = {xtj.ytpZt^). According to Pontryagin's minimum principle, X{tf) = 
■§^Zt\t=tf From this, it is straightforward to verify that (Ai(t/), A2(t/), X'i{tf)) = (0, 0, 1). 



Now let us consider another necessary condition A(t) 
the following relationships: 



A(t) 



(r(t),e(t),A(t),t) 
dr 



I Ai(t) ^ 




( Q -I \ 






1 -e 


\ Ht) ) 




^ e ; 




which leads to 



(35) 
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where {\i{t j) , \2{t f) , \^{t f)) = (0,0, 1). The corresponding solution to is 
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7I^sin^(t/ t) 




v 




J 




V 


lU^cosuitf t) + 1+^-2 


/ 



We obtain 

h{t)yt - ^2{t)zt = — ^[sincut + f sin cut/ + smio{tf - t)]. (37) 

It is easy to show that the quantity {\3{t)yt — X2{t)zt) occurring in fl32|) does not change 
its sign when tf E [0, :y^p] and t E [0, t/]. That is to say, e{t) = e is the optimal control 
when t e [0, ^]. Hence = zMt)) > Zt{i) = z\. 
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